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We study in general the time-evolution of correlation functions in a extended quantum system after 
the quench of a parameter in the hamiltonian. We show that correlation functions in d dimensions 
can be extracted using methods of boundary critical phenomena in d + 1 dimensions. For d — 1 
this allows to use the powerful tools of conformal field theory in the case of critical evolution. 
Several results are obtained in generic dimension in the gaussian (mean-field) approximation. These 
predictions are checked against the real-time evolution of some solvable models that allows also to 
understand which features are valid beyond the critical evolution. 

All our findings may be explained in terms of a picture generally valid, whereby quasiparticles, 
entangled over regions of the order of the correlation length in the initial state, then propagate 
with a finite speed through the system. Furthermore we show that the long-time results can be 
interpreted in terms of a generalized Gibbs ensemble. We discuss some open questions and possible 
future developments. 

I. INTRODUCTION 

Suppose that an extended quantum system in d dimensions (for example a quantum spin system), is prepared at 
time t = in a pure state \ip$) which is the ground state of some hamiltonian Hq (or, more generally, in a thermal 
state at a temperature less than the mass gap to the first excited state) . For times t > the system evolves unitarily 
according to the dynamics given by a different hamiltonian H, which may be related to Hq by varying a parameter 
such as an external field. This variation, or quench, is supposed to be carried out over a time scale much less than 
the inverse mass gap. How does the state \4>(t)) — e ~ lHt \4 l o) evolve? For a finite number of degrees of freedom, 
the system generically shows a periodic (or quasiperiodic) behavior with a period that typically increases when the 
number of degrees of freedom grows. This is the well-known phenomenon of quantum recurrence. However, in the 
thermodynamic limit this is no longer necessarily the case, and the natural question arises as to whether the system 
(or, rather, a macroscopically large subsystem) reaches a stationary state for very large times. To attack this question 
we consider the simpler question of how the correlation functions, expectation values of products of local observables, 
evolve and whether they reach constant values for large times. 

This problem has its own theoretical interest, being a first step towards the understanding of equilibration in 
quantum systems. On the same fundamental level we can also ask whether the approach can provide a new tool for 
the characterization of the collective excitations in strongly correlated systems. However, until recently it has been 
considered a largely academic question, because the time scales over which most condensed matter systems can evolve 
coherently without coupling to the local environment are far too short, and the effects of dissipation and noise are 
inescapable. Recent developments of experimental tools for studying the behavior of ultra-cold atoms have revised 
completely this negative attitude. In fact, thanks to the phenomenon of Feshbach resonance it is possible to tune the 
coupling of an interacting system to essentially any value and in extremely short times. In addition, the coupling to 
dissipative degrees of freedom can be much weaker than in ordinary condensed matter systems. This has allowed the 
observation, among other things, of the collapse and revival of a Bose-Einstein condensate [l[ and of the quenching 
of a spinor condensate [2j]. Furthermore, using highly anisotropic optical lattices, it has been possible to build and 
study essentially one-dimensional systems Q, and, of direct interest for this paper, the coherent non-equilibrium 
dynamics of these integrable models has been measured (Other interesting experiments are considered in Refs. 
[5(.) These striking experimental results have largely motivated the development of new numerical methods to study 
non-equilibrium dynamics, the most successful one being the time-dependent density matrix renormalization group 
(DMRG) @. 

On the purely theoretical side these kind of questions were first considered (as far as we are aware) in the seventies 
in the context of the quantum Ising-XY model in Refs. @, H[ (see also [1, EH El EH, EH EH for recent developments). 
However, in the very last few years, after the previously mentioned experimental progress, the study of quantum 
quenches has been pursued in a systematic manner and a large number of results is nowadays available. To quote only 
a few examples, the models considered include several realizations of one-dimensional (ID) Bose gases [TBI. Il6l. Il7l. fl8l. 
El, [H HI [H [H, H3, H, Luttinger liquids [H, [27j], coupled ID condensates [Hj], strongly correlated ID fermions 
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[29(, and mean field fermion condensates [3Ct l3l|. A closely related topic, that will not be considered here, concerns 
the formation of defects when crossing a critical point with changing the external parameter at a fixed rate, i.e. the 
so-called quantum Kibble-Zurek mechanism [32| . 

Most of these papers concern the exact or approximate solutions of very specific models. These are often strongly 
relevant, since in many cases they can be directly compared with existing experimental results or give very accurate 
predictions for future investigations. However, from the study of specific models it is difficult to draw general con- 
clusions on the physics of quantum quenches. There have been at least two notable exceptions. In Ref. [l5| it was 
conjectured that the asymptotic state at very large times can be described by a generalized Gibbs ensemble (for more 
details we refer the reader to section I VII where this conjecture will be explicitly discussed in relation to some of our 
findings). In Ref. [33[ we studied the time evolution after a quench in general, exploiting the path integral approach 
and the well-known mapping of the quantum problem to a classical one in d + 1 dimensions. The translational in- 
variance in the (imaginary) time direction is explicitly broken and thus the initial state plays the role of a boundary 
condition. This limits the range of applicability of field-theoretical methods, but when the hamiltonian H is at or 
close to a quantum critical point we could use the renormalization group (RG) theory of boundary critical behavior 
(see, e.g., [34|V From this point of view, particularly powerful analytic results are available for d = 1 because then 
the 1 + 1-dimensional problem is described asymptotically by a boundary conformal field theory (BCFT) [351 136| | . 

The aim of this paper is twofold: on the one hand we give a detailed description of the methods and the results 
reported briefly in Ref. 33], and on the other we generalize these results and give some new physical insights. Several 
results, in fact, appear here for the first time. The paper is organized as follows. In Sec. [II] we introduce the path- 
integral formalism that is applied in the following section HTT1 to one-dimensional critical systems by means of CFT. In 
Sec. lIVI we consider two simple ID models whose non-equilibrium dynamics is exactly solvable to check the correctness 
of the previous results and to extend our findings to gapped and lattice systems. In Sec. [V] we generalize the method 
to higher dimensions. In the last section [VI] we analyze all our results and explain most of the general findings in 
quantum quenches by means of simple physical arguments. We broadly discuss some open questions and possible 
future developments. 



II. PATH INTEGRAL FORMULATION AND SURFACE CRITIC ALITY 



Let us consider a lattice quantum theory in d space dimensions. The lattice spacing is a, and the lattice variables 
are labelled by a discrete vector variable r. Time is considered to be continuous. The dynamics of the theory is 
described by the hamiltonian H. Suppose we prepare this system in a state l^o) that is not an eigenstate of H and 
unitarily evolve it according to H . The expectation value of a local operator ©({r^}) at time t is 

(0(t, {rj)> = (Me im O{{Ti})e- iHt \^) . (1) 

We modify this time-dependent expectation value as 

(0(t, {r,})) = Z- 1 (^ |e^ t -^0({rJ) e -^ t -^|V.o) , (2) 

where we have included damping factors e~ eH in such a way as to make the path integral representation of the 
expectation value absolutely convergent. The normalization factor Z = (4>o\e~ 2eH \ipo) ensures that the expectation 
value of the identity is one. At the end of the calculation we shall set e to zero. 

Eq. ^ may be represented by an analytically continued path integral in imaginary time over the field variables 
</>(r, r), with initial and final values weighted by the matrix elements with \ipo): 

J [#(T,r)]^o|^(r 2 ,r))^(T 1 ,r)|^o)e-^ I ' MdT (3) 

where L[(f)]dT is the (euclidean) action. The operator O is inserted at r = 0, and the width of the slab is 2e. t\ 
and Ti should be considered as real numbers during the calculation, and only at the end should they be continued to 
their effective values ±e — it. In this way we have reduced the real-time non-equilibrium evolution of a d dimensional 
systems to the thermodynamics of a d + 1 field theory in a slab geometry with the initial state IV'o) pl&yiug the role 
of boundary condition at both the borders of the slab. The validity of our results relies on the technical assumption 
that the leading asymptotic behavior given by field theory, which applies to the Euclidean region (large imaginary 
times), may simply be analytically continued to find the behavior at large real time. 

Eq. ^ can be in principle used to describe the time evolution in any theory and for all possible initial conditions. 
Unfortunately, in confined geometries only a few field theories with very specific boundary conditions can be solved 
analytically in such a way to have results that can be continued from real to complex values. Some examples of these 
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will be given in Sec. [Vj (An interesting case concerns integrable massive boundary field theories as for example the 
recent application to the dynamics of coupled ID condensates in Ref. [I^].) However, the treatment greatly simplifies 
if a system is at or close to a (quantum) phase transition. In fact, in this case, we can use the powerful tools of 
Renormalization Group (RG) theory of boundary critical phenomena (see, e.g. [34|). For example, in the case of a 
scalar order parameter (corresponding to the Ising universality class), the slab geometry of above is usually described 
by the action [34[ 




where m 2 measures the distance from criticality, g > ensures the stability in the broken-symmetry phase, and c 
denotes the surface enhancement, i.e. the difference of the interactions on the surface with respect to the bulk. The 
value of c depends on the boundary conditions. For example, the choice c = +00 forces the field to vanish on the 
boundary, thus corresponding to Dirichlet boundary conditions. On the other hand c — —00 forces the field to diverge 
at the boundary. Finite values of c correspond to intermediate boundary conditions. According to the RG theory, 
the different boundary universality classes are characterized by the fixed point of c. A complete RG analysis (see, 
e.g. ,[13]) shows that the possible fixed points of c are c* = 00, 0, —00, with c* = ±00 stable and c* = unstable. Any 
other value of c flows under RG transformations to ±00. These universality classes are called ordinary (c* = 00), 
special (c* = 0), and extraordinary (c* = —00) [34[. 

Thus, for the purpose of extracting the asymptotic behavior, as long as \tpo) is translationally invariant, we may 
replace it by the appropriate RG-invariant boundary state \i/jq} to which it flows. The difference may be taken into 
account, to leading order, by assuming that the RG-invariant boundary conditions are not imposed at r = t% and 
Ti but at t = T\ — To and r = t-i + tq. In the language of boundary critical behavior, tq is called the extrapolation 
length [34[ • It characterizes the RG distance of the actual boundary state from the RG-invariant one. It is always 
necessary because scale-invariant boundary states are not in fact normalizable[36j]. It is expected to be of the order of 
the typical time-scale of the dynamics near the ground state of Ho, that is the inverse gap m^" . (This can be checked 
explicitly for a free field theory, see later). The effect of introducing To is simply to replace e by e + To. The limit 
e — > 0+ can now safely be taken, so the width of the slab is then taken to be 2to. For simplicity in the calculations, in 
the following we will consider the equivalent slab geometry between t = and r — 2tq with the operator O inserted 
at t = To + it. This is illustrated in the left part of Fig. [TJ In fact, in this geometry we can consider products of 
operators at different times tj , by analytically continuing their labels to To + itj . 

Let us discuss to what initial conditions the fixed points correspond. In the case of free boundary conditions, when 
the order parameter is unconstrained, the continuum limit forces it to vanish there, so it corresponds to the ordinary 
transition (c = 00). For a lattice model, this corresponds to a completely disordered initial state (e.g. for the Ising 
chain in a transverse field corresponds to infinite field). In contrast for non- vanishing fixed boundary conditions (e.g. 
+ or — for Ising- like systems), the continuum limit makes the order parameter to diverge at the boundary, thus 
corresponding to the extraordinary transition (c = —00). 

III. ONE SPACE DIMENSION AND CONFORMAL FIELD THEORY 

In this section we specialize the methods just introduced to the case when H is at a quantum critical point whose 
long-distance behavior is given by a 1+1-dimensional CFT, with dispersion relation f2/- = v\k\. We set v = 1 in the 




FIG. 1: Left: Space-imaginary time region in ([2]). Imwi = r, that will be analytically continued to r — > to + it. Right: 
Conformal mapping of the left geometry to the upper half-plane (c.f. Eq. ©). Note that argZi = 9 = 7tt/2tq. 
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following. RG-invariant boundary conditions then correspond to conformally invariant boundary states. In this case 
the correlation functions are accessible through the powerful tools of boundary CFT. 

The main property we will repeatedly use in the following is the relation of correlation functions of (primary) 
operators among two geometries connected by a conformal transformation. For example, the slab geometry of above 
is just a two-dimensional strip whose points are labelled by a complex number w = r + it with < Imw < 2tq. The 
strip can be obtained from the upper half-plane (UHP) Imz > by the conformal mapping 

w{z) = — - logz , (5) 

with the images of points at the same imaginary time on the strip lying along argz^ = 6 = ttt/2to. In the case where 
O is a product of local primary scalar operators <i>i(wi), the expectation value in the strip is related to the one in the 
UHP by the standard transformation 

<n $ iM>strip =Y[\w'(Zi)\- Xi (Q^O^UHP , (6) 

i i i 

where Xi is the bulk scaling dimension of Note that the (eventual) expectation values of the $i in the ground state 
of H are supposed to have been subtracted off. The asymptotic real time dependence is obtained via the analytic 
continuation t — * tq + it, and taking the limit t, 3> tq. 

In the following subsections we apply these methods to some specific cases. 



A. The one-point function 

In the UHP, the one-point function of a scalar primary field with bulk scaling dimension x is (3>(.z))uhp = 
Af [2lm(z)]~ x , as a simple consequence of scaling invariance. The normalization factor Af is a non-universal ampli- 
tude. In CFT the normalizations are chosen in such a way that the bulk two-point functions have unit amplitude 
(i.e. ( < I > (2i) ( I ) (22))buik = \%2 — zi\~ 2x ). This choice fixes unambiguously the amplitude Af that turns out to depend 
both on the considered field $ and on the boundary condition on the real axis b. It vanishes if $ corresponds to an 
operator whose expectation value in \ipo) vanishes, and thus ($(£)) = 0, for all times. 

When the primary field is not vanishing on the boundary, performing the conformal mapping J5| we obtain 



($H) 



strip 



\ w '(z)\-*m Z (w))) 



UHP 



A b 



1 



4r sin(7rr/(2T )) 



(7) 



that continued to real time r = r + it gives 



(*(*)>= A 



l 



4t cosh(7rf/(2T )) 



(— 

\2n 



D — x7rt/2ro 



(8) 



Thus the order parameter (and any other observables described by a primary field) decays exponentially in time to 
the ground-state value, with a non- universal relaxation time t®, = 2to/xqtt. The ratio of the relaxation times of two 
different observables equals the inverse of the ratio of their scaling dimensions and it is then universal. 

The normalization factor Af is known for the simplest boundary universality classes (37| . In the case of $ being 
the order parameter and the boundary condition is fixed (ipo(z) — °°) Af is 1 for the free boson and 2 1 / 4 for the 
Ising model. 

An important exception to this law is the local energy density (or any piece thereof). This corresponds to the tt 
component of the energy-momentum tensor T^ v . In CFT this is not a primary operator. Indeed, if it is normalized 
so that (T^tjhp = 0, in the strip (38| (T ti (r, r)) = 7rc/24(2r ) 2 (where c is the central charge of the CFT) so that 
it does not decay in time. Of course this is to be expected since the dynamics conserves energy. A similar feature is 
expected to hold for other local densities corresponding to globally conserved quantities which commute with H , for 
example the total spin in isotropic models. 



B. The two-point function 



In the case of the one-point function, scaling invariance was enough to fix the functional dependence on the position 
in the UHP. However, the form of the two-point function depends explicitly on the boundary universality class and 
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on the operator considered. In the following subsections we will consider the equal-time correlation function for the 
order parameter in the gaussian and in the Ising universality classes that are easily treated in full generality. At the 
quantum level they describe (among the other things) a chain of harmonic oscillators (explicitly considered in Sec. 
IIV A[) and the Ising model in a transverse field (whose real time evolution has been considered in Refs. @) H, H El and 
is briefly reviewed in Sec. IIV Bp . Finally we will discuss the general form of the two-point function for asymptotically 
large time and distance, that can be obtained from general CFT arguments. 

1. The gaussian model 

The content of this subsection has been already reported in Ref . [39( , during the study of the time evolution of the 
entanglement entropy, that transforms like the two-point function of a primary field in a boundary gaussian theory 
[40| . We report it here for sake of completeness. 

For a free boson the two-point function in the UHP is [3{| 



(*(*i)*(*a)>UHP = , (9) 

\Zl2Z 12 Z 11 Z 22 / 

with Zij — \zi — Zj \ and zj, — "z/T. Note that <& is not the gaussian field 0(z), but its exponential $(z) — e l9 ^ (see Sec. 

Under the conformal mapping ([5]) we obtain the two-point function on the strip at imaginary time r, at distance r 
apart 

(S(r,T)3(0,T)) Btrip = \w'(z 1 )\- x \w , (z 2 )\- x mz 1 (w))mz 2 (w))) mP = 

cosh(7rr/2ro) — cos(ttt/tq) 



2t / 8 sinh (irr /4r ) sin 2 (ttt /2t ) 



(10) 



that continued to real time r = tq + it gives 



(*(r,*)$(0,t)) 



7r \ cosh(7rr/2ro) + cosh(7rt/ro) 



2t J 8sinh 2 (7rr/4r ) cosh 2 (7ri/2r ) 
In the case where t/tq and t/ro are large this simplifies to 



(11) 



$ (r,t $ (0,t )) = (n 2t ) 2x [ /n w cx < ,„ , . (12) 

\ V . ; V ) 7/ WO; ^ e 7rr/2r . e 7tt/r J |g-x7rr/2r fort>r/2 V ' 

i.e. the two point function at fixed r decays exponentially in time up to t* — r/2 and then saturates to a value that 
depends exponentially on the separation. 

In the case of fixed initial condition, with one-point function given by Eq. l[8"]), the connected correlation function is 

(*(r,t)*(0,t)U s = (Hr,t)mt)) - (*(0,t)) 2 « r_ W2T0 _ e ^ t/T0 r J/ 2 \ (13) 

i.e. correlations start developing at t* = r/2 and, being t 3> To, at t* the connected two-point function almost 
immediately jumps to its asymptotic value. In the case of disordered initial conditions (ipo{r) — 0), connected and 
full correlation functions are equal. 

2. The Ising universality class 
For the Ising model the two-point function in the UHP is (3f| 

($(z 1 )$(2 2 )) UHP - ( * 15 * 21 ) ^ F(rj) , (14) 
\z\ 2 z 12 z Y1 z 22 J 
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where F(rj) is given by 



F( V ) = 



Vj + i] 1 / 2 ±y/l- rf/ 2 

71 



and rj is the four-point ratio 



V 



z \l z 22 



(15) 



(16) 



The sign ± depends on the boundary conditions. + corresponds to fixed boundary conditions and and — to disordered 
ones. 

The only difference with respect to the gaussian case is that we have also to map F(rj) according to the conformal 
transformation ([5]) . After simple algebra we have 



2sin 2 (7rr/2r ) 



cosh(7rr/2ro) — cos(7tt/to) 



and so 



($(r,r)$(0,r)) 



strip 



7r \ cosh(irr / 2tq) — cos(itt / tq) 
2t J 8sinh 2 (7rr/4To)sin 2 (7rr/2To) 



1 1/8 



1 

71 



y/2 sin(7rr/2To) 



\/cosh(7rr/2ro) — cos(7tt/to) 
Analytically continuing to real time t = tq + it we obtain 



±4/1 



y/2 sin(7rr/2ro) 



-\/cosh(7rr/2To) — cos(7rr/ro) 



($M)$(0,t)> = 



cosh(7rr/2r ) + cosh(7rf/ro) 
2t J 8sinh 2 (7rr/4T ) cosh 2 (7rt/2r ) 



1/8 



1 

71 



\/2cosh(7rt/2ro) 



-\/cosh(7rr/2To) + cosh(7rf/To) 
that for t/tq and t/ro much larger than 1 simplifies to 

i/ 4 1 / e 7rr / 2T ° + e 7r */ T ° \ 
(7r/2r ) -j= £7rr/2T0 . eJrt/T0 ) 



±4/1 



v / 2cosh(7rt/2r ) 



■y/cosh(7rr/2ro) + cosh(7ri/r ) 



e 7rt/2r 



^/ e irr/2r ZjZ g7rt/T 



±W1- 



^/ e 7Tr/2T _|_ g7Tt/T 



(17) 



(18) 



(19) 



(20) 



Note that the exponential terms in the square root are always -C 1. Thus for fixed boundary condition we get the 
free boson result Eq. (fT2"]) with 2 = 1/8. For the connected part we need to subtract ($(0,i)) 2 given by Eq. ([5]) with 
Af = 2 1 / 4 . We finally obtain 



(*(r,i)3(0,i)) t 



($(r,0$(0,i))-($(0,t)) 2 cx 



for£<r/2 : 

e -7rr/16r _ g-Trt/STo for t>r f 2 . 



(*(r,t)*(0,t)) cx 



r/2r 1 p Tvt/Ta 



1/8 



O 7rt/2T 



e 7rr/2r . gTrt/ro / ^/ e nr/2r ZjZ e 7rt/T 



e -7r(r-3/2i)/4r fori < r/2, 

fort > r/2 , 



-7rr/16ro 



(21) 



Thus, also for the Ising model with fixed boundary conditions, connected correlations start developing at t = i* = r/2. 
In the case of disordered initial condition, we have 



(22) 



resulting in an exponential space dependence even for t < r/2 (clearly in this case the connected correlation function 
equals the full one). 
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3. The general two point- function 



From the results reported for the gaussian and Ising models, it is now relatively simple to understand the general 
properties of the time dependence of the two-point function in the very general case. The two-point function in the 
half-plane has the general form [35| 

<$(z 1 )<i>(z 2 )) U HP = ( Zl * Z& ) X F(ri) , (23) 

where the function F(rf) depends explicitly on the considered model. Under the conformal map to the strip we know 
that the first part of Eq. ([23| transforms according to Eq. (jTOj) . Thus we need only to map F(r]) that, in the 
general case, is an unknown function. However, during the study of the Ising model we showed that the analytical 
continuation of r) for t, r ^ To is 

fj r^j - (24) 

' e nr/2r n _|_ gTrt/r ' \ ' 

Thus for t < r/2 we have r\ ~ e 7r (*- r / 2 )/ T o <g 1 anc j i n the opposite case t > r/2 we have r\ ~ 1. As a consequence 
to have the asymptotic behavior of the two-point function we only need to know the behavior close to r\ ~ (i.e. 
the behavior close to the surface) and for 77 ~ 1 (i.e. deep in the bulk). Fortunately they are both exactly known. 
Indeed when 77 ~ 1 the two points are deep in the bulk, meaning F(l) = 1. Instead for rj -C 1, from the short-distance 
expansion, we have 

F( v )c(Af)^, (25) 

where Xb is the boundary scaling dimension of the leading boundary operator to which $ couples and Af is the 
bulk-boundary operator product expansion coefficient that equals the one introduced in Eq. {8| [see e.g. Ref. [37|]. 
All the previous observations and the explicit calculations of the previous sections lead for t > r/2 to 

($(r,t)$(0,<)) oc e -W2r ; (26) 

while for t < r/2 we get 

($(r, t)*(0, t)) cx {Af f e - x ^ T ° x e **»(.*-r/2)/n . (27) 

Note that if ($) 7^ 0, x;, = and the last factor is absent. The leading term is then just ( < I > ) 2 . Thus the leading 
term in the connected two-point function vanishes for t < r/2, and its first non- vanishing contributions is given by 
subleading terms either in F or in the bulk-boundary short-distance expansion. 

It is very interesting that we only have to know the behavior as r\ — > and 1 to get the results we need for large r 
and t. However, we stress that only a complete calculation (as those performed in the preceding sections) gives the 
full analytic structure of the CFT result needed to justify the analytical continuation from imaginary to real time. 
Moreover, the behavior within a distance 0(tq) of the horizon r = 2t depends on the detailed form of F. 



C. Correlations functions at different times 



Let us consider the case of the two-point function calculated at different real times ($(r, t)$(0, s)). This is again 
obtained by mapping the imaginary time strip to the UHP, but in this case the two points are w± = r + ir\ and 
W2 = + 7>2, that, at the end of the calculation, must be analytically continued to t± = tq + it and Ti = tq + is. See 
Fig. [2] for a pictorial representation of the space-time domain in the strip and the resulting mapping to the UHP. 

Let us start the discussion with the free boson. The distances appearing in Eq. ^ are 

zf 2 = 1 + e 7rr/To - 2e nr/2T0 cos(0i - 6 2 ) , zn = 2 sinfli , 

zl- 2 = 1 + e nr / ro - 2e 7rr / 2ro 008(61! +6 2 ), z 22 = 2e 7!r ' 2To sin6» 2 , (28) 

where Qi — 7TTj/2To. Thus the correlation function on the strip is 

($(r,Ti)$(0,T 2 )) strip = \w'(z 1 )r x \w'(z 2 )\- x mz 1 (w))mz 2 (wj)) VHP = (29) 

/ tt \ 2 cosh(-7rr/2T ) - cos(7t(ti + t 2 )/2t ) 

\2t J 4sin(7rri/2T ) sin(7rT 2 /2T )(cosh(7rr/2T ) - cos(7r(r! - r 2 )/2r )) 
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FIG. 2: Left: Space-time region for the correlation functions at different times. Right: Conformal mapping to the upper 
half-plane. 



that for Ti = T2 reduces to Eq. (fTO)) as it should. Continuing to real times and considering r, t, s, \t — s\ ^> tq we 
obtain 



,-x7r(i+s)/4ro for r > t + S . 



STu S | -S e ' W4T ° fort- S <r<i + S , (30) 

e . v ; feT0+gT0 i n | e _ xff | t _ 5 | /4ro forr < It - s| . 



Following the line sketched in the previous subsection, it is easy to generalize this result to the most general CFT. 
In the case of a theory with fixed initial conditions (i.e., ($) ^ 0) the asymptotic result is the same as before, with 
only the crossover points being affected by the precise expression for F(rj). Instead, in the case where (<&) = 0, the 
first case gains an additional factor e~ 7ra5i> ( t+s_r )/ 4r ° _ 

Note that the autocorrelation function (i.e. r = 0) has only an exponential dependence on the time separation t — s 
and does not exhibit aging in this regime. 



D. Evolution with boundaries 



We now consider the case of time evolution of a half-chain with some boundary condition at r = 0. For simplicity 
we assume that the (conformal) boundary condition is of the same kind of the initial boundary condition (for example 
we fix all the spins at t = and at the boundary r = to point in the same direction). The space-time region we 
have to consider is depicted in Fig. [3] If different initial and boundary were considered, one needs to insert boundary 
conditions changing operators at the corners of the figure. 

The w plane is mapped into the UHP by 



. TTW 

z{w) = sm — , 

ZTq 

with the corners at ±tq mapped to ±1. The mapping of w\ is 

z\ = z(w\) = z(—tq + t\ + ir) = — cos(7tti/2to) cosh(7rr/2ro) + i sin(7rri/2ro) sinh(7rr/2To) 
In the z plane the 1-point function is 

($0?i))uhp oc \ln\ Zl \- x -> [sin(7rTi/2To)sinh(7rr/2T )] _:l; , 



and 



2to 

7T 



1-Z 2 



OC 



\W'( Z1 )\ 2 : 

Thus on the strip we have 

($(wi)) str ip = |w'(zi)r ;r ($(w(2:i)))uHP oc 



cosh(7rr/ro) — cos(7tti/to) 



sin 2 (7rri/2ro) sinh 2 (7rr/2ro) 



cosh(7rr/ro) — cos(7tti/to) 



r/2 



that continued to real time t\ = it is 



(*(t,r)) 



cosh(7rf/ro) + cosh(7rr/ro) 



cosh("7rf/2T ) 2 sinh (7rr/2ro) 



x/2 



(31) 
(32) 
(33) 

(34) 
(35) 
(36) 
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FIG. 3: Left: Space-time region for the one-point function in a boundary (at r = 0) geometry. Note w = r + ir. Right: 
Conformal mapping to the upper-half plane, c.f. Eq. (|31[) . 



and for t, r 3> To simplifies to 



($(t,r) oc 



r/r 



x/2 



^ — ■nxr/2To 



for £ < r , 
for t > r . 



(37) 



Note that in this case the characteristic time is t* = r and not r/2. This explains also why the entanglement entropy 
of a semi-infinite chain with free boundary condition at x = 0, has characteristic time t* = r as firstly noted in Ref. 



E. Discussion and interpretation of the CFT results 



All the correlation functions calculated so far display two very general features: first there is a sharp horizon (or 
light-cone) effect at t = t* = r/2 (or r) resulting in a behavior before and after t* completely different; second the 
asymptotic long-time correlation functions are the same as those at finite temperature /3 e ff = 4to. The light-cone 
effect is a very general phenomenon and and will be discussed in section WT\ In the following sections we will point 
out that also the "effective temperature" is a general phenomenon, but it has some specific CFT features that are 
worthy of comment. 

In fact, it is very easy to understand the technical reason why we find an effective temperature despite the fact that 
we are studying a pure state at T = 0. The finite temperature correlations can be calculated by studying the field 
theory on a cylinder of circumference (3 = 1/T. In CFT a cylinder can be obtained by mapping the complex plane 
with the logarithmic transformation f3/{2it) \ogz. The form for the two-point function in the slab depends in general 
on the function F{rj) -cf. Eq. (fT5)) - but when we analytically continue and take the limit of large real time, we find 
that effectively the points are far from the boundary, i.e. at r] = 1. Thus we get the same result as we would get if we 
conformally transformed from the full plane to a cylinder, and from Eq. ([5]) the effective temperature is /3 c ff = 4to. 
A similar argument can be worked out for the multi-point functions as well. 

An addtional comment concerns what we expect for correlation functions of general operators, not only primary. 
We have clearly seen that the local energy density does not relax, of course, and this is consistent with its not being 
primary. But, we also know that at finite temperature (Ttt)p = irc/6j3 2 , that is perfectly compatible with the previous 
result with /3 e s = 4to. Furthermore, for large real times, the two-point function of T tt — (T) tt does behave as though 
it was at finite temperature. This means, in particular, that the energy fluctuations in a large but finite volume (the 
specific heat) are the same as those at finite temperature. Thus one is tempted to extend this finite temperature 
interpretation to non-primary operators, and indeed it is the case. In fact, in the argument of above for the equivalence 
of the long-time correlations and finite T, there are essentially three steps. First, we need to write down the form of 
the correlation function in the half-plane, but this only depends on special conformal transformation and so is valid 
for any quasi-primary operator like Second we have to transform from the slab to the half-plane: for non-primary 
operators this can have some anomalous term. Finally we need to compare the large limit of this to what one would 
get transforming directly from the cylinder to the full plane. However, the two conformal transformations we are 
comparing are both logarithmic (in one case 2ro/-7rlogz, in the other case (3/2nlogz), so the anomalous terms should 
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be the same. Thus the two correlations functions are the same. This argument works for all quasi-primary operators. 
Then, since we can get all the non-quasiprimaries by considering successive operator product expansions with the 
stress tensor, it also works for all operators. 



IV. EXACT REAL-TIME DYNAMICS IN SIMPLE INTEGRABLE CHAINS 

The results of the previous section rely on the technical assumption that the leading asymptotic behavior given 
by CFT, which applies to the euclidean region (large imaginary times), may simply be analytically continued to 
find the behavior at large real time. While such procedures have been shown to give the correct behavior for the 
time-dependent correlations in equilibrium, it is important to check them in specific solvable cases for non-equilibrium 
evolution. 

Thus, as a complement to the CFT calculations, in this section we consider the real-time evolution of two simple 
analytically tractable models. We solve the dynamics of a chain of coupled harmonic oscillators and we review and 
re-analyze some known results for the Ising-XY chain in a transverse magnetic field. Beyond providing examples of 
the CFT results (with central charge c = 1 and 1/2 respectively) in the critical case, these models allow us to take 
into account the effects of a finite mass-gap and of the lattice. 



A. The chain of harmonic oscillators 



The simplest model with an exactly solvable non-equilibrium dynamics is surely a chain of coupled harmonic 
oscillators with hamiltonian 



H = \ E + + E^+, - 0,) 2 ] • (38) 



We introduce a coupling more general than simple nearest-neighbor hopping so as to allow for a general dispersion 
relation below. For simplicity we also assume periodic boundary conditions. ip n and 7r„ are the position and the 
momentum operators of the n-th oscillator, with equal time commuting relations 

[<Pm, K n ] = iSnrn [<Pn, Pm] = fan, Km] = • (39) 

The hamiltonian can be written in diagonal form H(m) = J2k ^kA\Ak with modes 

Ak = —==(flkPk + i^k) , (40) 

ttl = TO 2 + 2^^(l-cos(27rfcj/7V)) . (42) 

3 

Note that we use the same symbols for the operators and their Fourier transforms (<pk — 1/VN J2n=o e 27rlkn ^ N <f n 
and analogously for 7Tfe). 

We consider the scenario in which the system is prepared in a state \tpo), that is ground-state of H(m ), and at the 
time t = the mass is quenched to a different value m ^ m . We use the notation Sl k for the dispersion relation for 
t < and the fife for the one for t > 0. 

Since (V'olVnlV'o) = 0; the expectation value of the field ip n vanishes at any time. This example in fact corresponds 
to the quench from the disordered phase in the language of the previous section. Thus we concentrate our attention 
on the two-point function 

W)\<Pn<P0W)) = W-Ol^n (*)¥>? WW , (43) 

where we introduced the operator in the Heisenberg picture <Pn(i), whose time evolution is given by 

N-l 



k=0 



V" (t) = E v m~ k ( ei{Pkn ~ nkt)A k + e -^"-°^4) , (44) 
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where p k — 2nk /N. Accordingly, the product of the two fields is 



■ 



(45) 



In order to have the time dependent two-point function we need the expectation values of the bilinear combinations 
of A's on the initial state, that is annihilated by the A$ k . Thus it is enough to write A's as functions of Ap's, i.e. 



A k = i 



.4 



4 



A 0k 



n k 

^Ofc 

fife 

^Ofe 



1 0-fc 



= CfeAofe + rffeAj_ fe , 



= CkAl k + d k A - k , 



leading to (we understand (•) = (V>o| ■ \ipo)) 

(A k A k >) = ((c k A ok + d k Al_ k )(c k >A ok > + d k >A\_ k ,)) = c k d k > {A ok A\_ k ,) = c k d k 8 k _ k , 

{A k A\,) = ((c k A ok + d k A\_ k ){c k r A\ k , + d k >A Q _ k ,)) = c k c k > {A ok A\ k ,) = c 2 k 5 ktk > , 

(A\A k >) = {{c k A\ k + d k A - k )(c k <A okl + d k 'Al_ k ,)) = d k d k , (A - k Al_ k ,) = d 2 k 5 k ^ , 

( A t A l>) = (i c k A lk + d kA - k )(c k/ Al k , + dfc/Ao-kO) = d kCk'(A - k Al k ,) = c k d k S k ^ k , 
Finally we arrive at 



(<p% (t)pf (t)) = 4~ E 7T \c k d k e^ n ' 2n ^ + c£e** w + die-™ + c k d k e-^ r 
iv a ' s si. . 

which, in the thermodynamic limit (iV — > oo), may be written as 



(<p?(t)<p»(t)) - (^(0)^(0)) 



^ cifcr (ng fc -ng)(i-cos(2n fc t)) 



BZ 



(2tt) 



(46) 
(47) 

(48) 
(49) 
(50) 
(51) 

(52) 
(53) 



where the integral is on th first Brillouin zone |fc| < it /a. Note that for t = and for m = mo this two-point function 
reduces to the static one, as it should. This result can also be found by integrating the Heisenberg equations of 
motion for each mode. For future reference it is also useful to write down explicitly the Fourier transform known as 
momentum distribution function 



p(k) 



(nl k + nt)-(n 2 0k -nl)cos(2n k t) 



'fei'Ofe 



(54) 



Note that when considering correlation functions in momentum space, a long-time limit does not exist and we need 
to take the time average, in contrast to what happens in real space. 



The continuum limit 



In Eq. (|52p everything is completely general and applies to any chain with finite lattice spacing. Let us know discuss 
the continuum limit that is achieved by sending N — > oo in such a way that (1/N) J2 k ~ * /_ dp/(2ir), p n — > p, and 
<>I • fl? m 2 ■ ,r. 



In this limit the correlation function becomes (Am 2 = m\ — m 2 ) 



dp 
— ( 

2tt 



ipr 



-Am 2 cos(2y/p 2 + m 2 t) + m 2 + m 2 , + 2p 2 
(to 2 + P 2 )y / nT'Q + p 2 



(55) 



Since a closed form for such integral is quite difficult to write down in the most general case, we will only consider 
some particular cases. 




0.2 0.4 0.6 0.8 1 2 4 6 8 

t t 

FIG. 4: Left: G(r,t)/mo given by Eq. I|60[l as function of t, at fixed r = 1. Three different values of mo = 10,3, 1 (from the 
bottom to the top) are shown. Inset: Lattice effects showing the cos4t oscillations on top of the continuum result. Right: 
G(r,t)/mo given by the numerical integral of Eq. (|62|l as function of t, at fixed m, r = 1. It is compared with the asymptotic 
behavior for < 2t — r <C m _1 and for i>r. 



Let us first consider the conformal evolution (to = 0) from an initial state with a very large mass too — > oo. This 
should reproduce the CFT result previous section for all the time t, since the correlations of the initial state shrink 
to a point. We have 

^ f 00 dp lvr l-cos{2pt) f° fori<r/2, 

^"^y^^— ? L ^- mo (i(2 t -r) fa*>r/2. (56) 

To compare such result with the conformal result given by Eq. (|12p. we have to keep in mind that the primary field is 
not the gaussian one ip H (r,t), but its imaginary exponential. Thus we need {e lqipH ( r ' 1 ) e~ lqipH (°>*)), with an arbitrary 
q. Despite of the apparent complexity of such correlator, it is very simple to obtain it using the standard property of 
gaussian integrals 



{e iqip l <(r,t) e -iq V H '(0,i)) = g -g 2 {(<p H (r,t)- v H (0,t)) 2 )/2 = & q 2 {G{r,t)-G(0,t)) ^ 



leading to 



, i,/(r,() (o,t)v _ I e 9 " mot fOT * < r/2 , , , 



that is exactly the same of Eq. (|12[) with x$ oc q 2 and ro oc to,q , confirming that To is just proportional to the 
correlation in the initial state, as its interpretation in terms of the extrapolation length suggests. 

The case to = and too finite corresponds to a conformal evolution from a generic state with correlations propor- 
tional to TOq 1 . Thus we expect the CFT result Eq. fTS]) to be true for asymptotic large times and separations. 

From Eq. (|55p. the two-point function of the gaussian field is 

GM) = r dp ripr - m lcos(2pt)+rn* + 2 P > 
J-oo 27r P 2 y/ml +p 2 



that can be written as 



on A ok ( U« ^ f(m (r-2t)) + f(m (r + 2t)) 

2nG(r, t) = 2K (m r) + f{m r) , (60) 



where Ko(y) is the modified Bessel function and 



1^21 (V 2 3/2 



f(y) = i + ~G? 3 r T 



1 1/2 



(61) 
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with Gjl tne Meijer G-function (see e.g. [42^). Note that f(x) is characterized by f"(x) = -K (x) and /(0) = /'(0) = 
0. In the limit too — * oo, it is easy to show that G(r,t) reduces to the previous result. In the left panel of Fig. 31, we 
report G(r, t) as function of t at fixed r — 1 for mo = 10, 3, 1. It is evident that a finite Too results in smoothing the 
curve close to t = r/2 and giving an offset in zero. Both the effects are more pronounced as too decreases. For large 
t, independently on mo ^ 0, we have G(r,t) = t + 0(t°), confirming that the CFT result is correct for asymptotic 
large times. 

The case with arbitrary m and mo is quite cumbersome to be worked out analytically and not really illuminating. 
For this reason we concentrate here on the massive evolution from a state with too — > oo. In this case the correlation 
function reads 



,. i d P i P r 1 _ cos(2yy + m 2 t) 
G(r,t) =m I —e p , 2 , ■ (62) 



2n (to 2 + p 

In the limit t — ► oo the cosine term averages to zero, giving 



-mr 



/°° dp rriQ g 
7^7 — o tT ~ « — ■ (63) 
.oo 2tt (m 2 +p 2 ) m 2 y ; 

To understand the time dependence, let us first note that, despite the presence of a square-root, the integrand of 
Eq. (|62[) is analytic in p, since the square-root is the argument of the even cosine function. Thus we can make the 
integral in the complex plane and use the Cauchy theorem. As long as r > 2t, the behavior for p — > ioo is dominated 
by e tpr and we can safely close the contour path in the upper half-plane where the residue at p — im is zero. As a 
consequence G(r,t < r/2) = 0. 

For t > r/2, the integral is more difficult, since we can not close the contour in the upper half-plane, because the 
cosine is "larger" than e lpr for p — > ioo. However, the approach to the asymptotic value is easily worked out. In fact, 
for t>r, the term e lpr in the integral (|62p can be approximated with 1, since it is slowly oscillating. Thus we have 



s s f dp cos(2 J p 2 + m 2 t) , , . 

G(r, t > r) - G(r, t = co ~ -m / -f \\ p T ' = -m f m (t) , 64 

where f m (t) = l/(2m) - t 1 F 2 (l/2; 1, 3/2; ~m 2 t 2 ), that satisfies f m (t) = -J (2mt) and / m (0) = l/(2m). Note that 
f m (t) is just the result for r = 0. In the complementary region < 2t — r -C to^ 1 , the integral in Eq. (|62|) is dominated 
by the modes with p>m, and so it can be described by the conformal result. 

In the right panel of Fig. 31 we plot the time-dependence of G(r, t) at fixed r as obtained by numerically integrating 
Eq. (|6"2"|) for to = r = 1. The plot shows that the conformal result t — r/2 describes the behavior close to t = r/2, 
while for larger times the asymptotic expression gives an excellent approximation. 

We finally note that G(r,t) at fixed t displays spatial oscillations (for r < 2t, else it vanishes), as it can be simply 
realized by a stationary phase argument. 



2. Lattice effects 



Another advantage of this simple model is that the effects of the lattice can be easily understood. In fact, by 
a stationary phase argument, the dominant contribution to (|53[) in the limits of large t and r comes from where, 
the group velocity Vk = — r/2t, independently of the explicit form of Consequently, the two-point function 
of a gapless model with dispersion relation fi^ = 2sinfe/2 differs from the continuum limit previously derived for 
t > r/2v m , where v m is the maximum group velocity. In particular G{r,t) receives a contribution from the slowest 
mode (k = ir with d t = 0) whose effect is to add to Eq. (|5"6")l fast oscillations going as cos^Sl,^) = cos4t. The 
resulting G(r, t) is plotted in the inset of Fig. 31 m the gapless case. 

Thus, lattice effects play a more important role in the cases where the asymptotic result vanishes. A typical example 
is the energy density that in this model is proportional to (<p r +i ~ Vr) 2 ■ The continuum calculation would just give 
zero (the mean-energy is conserved), but the approach to this value is governed by lattice effects and is dominated by 
the smallest group velocity that comes from the zone boundary at \k\ — tt. In fact, just with a trivial calculation we 
have 

»/ ,o, 9, , f dk , , . (1 — cos 2f2fct) 

(foVfi - Vr?) = 2«p§)-< m > 1 »= / — (l-cosfc) V • - 



2tt x ' 2sin 2 fc/2 

/ — (1 - cos2r2 fc f) oc J (4t) - t" 1/2 cos4i , (65) 
J 2n 
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that displays the typical cos4i from the zone boundary at |fe| = tt. 

The lattice dispersion relations are sensitive to the microscopical details of the model and consequently quantities 
like the energy-density show a dependence on this. For example, a lattice massless fermion has dispersion relation 
f2fc = 2sinfc/2 and energy density given by sin 2 k. In this case the time evolution of the energy density is 

* s . n2fc (l_ccj2 M 

2tt 2sin 2 fc/2 V ' 

that for large times goes like i -3 / 2 cos4t, resulting in a different power law compared to before. 

If the system is quenched to a gapped H with fi| = m 2 + 2(1 — cosfc) the maximum group velocity corresponds 
to a non-zero wave number. This gives rise to spatial oscillations in the correlation function (this is true also in the 
continuum limit as already discussed in the previous subsection). 



B. The Ising-XY chain in a transverse magnetic held 



The most studied one dimensional quantum spin model is the so-called XY chain in a transverse field, defined by 
the hamiltonian 



H 



N 

£ 



(! + 7) x 



°3 



(1-7) 



v y 



ha. 



(67) 



where a^' y ' z are the Pauli matrices, 7 is called anisotropy parameter, and h is the applied external transverse held. 
It is well known that for any 7 7^ the model undergoes a phase transition at h = 1, that is in the universality class 
of the Ising model (defined by 7 = 1). For simplicity we will just consider the Ising case, other values of 7 7^ being 
equivalent. 

We consider the non-equilibrium unitary dynamics that follows from a quench of the magnetic field at t — from 
ho to hi 7^ ho. Earlier works on this subject by Barouch et al. date back to seventies @, H| . Exploiting the mapping 
of this model onto a free fermion, the time evolution of the transverse magnetization (that is not the order parameter) 
was obtained exactly 0: 



/ \ 1 
m z [t) = — 



dk 



(ho — hi) sin 2 k cos(2eit) — (cos k — /ii)[(cos k — ft,o)(cos k — hi) + sin 2 k] 



eoe 2 



(68) 



where ej = e(/ij) with e(h) = \J sin 2 k + (h — cos k) 2 . In Ref. 0] it was shown that (for non-exceptional parameters) 

the approach to the asymptotic value for t — > 00 is of the form i -3 / 2 cos4£. This is simply shown in the case of ho = 00 
and hi = 1, when the integral simplifies and we obtain (see also [9() 



m z (t) = - + - 

Z 7T 



dk 



sin 2 k cos(2eit) 



Ji(4t) 



(69) 



,0 4sin 2 fc/2 2 At 

where J\ is the Bessel function of the first kind whose asymptotic expansion for large argument is J\(x) ~ 
y/2/irx cos(x + 7r/4). 

To understand this result we should keep in mind that the transverse magnetization is not the order parameter: 
it corresponds to the product of two "disorder parameters" at neighbor sites on the dual lattice. Thus it must have 
the symmetry of an energy operator, for which CFT just predicts a constant asymptotic result. The power law term 
t -3 / 2 cos(4£) is just the lattice correction to the asymptotic. This correction clearly shows the fermionic nature of the 
model on the basis of what discussed in the previous section. 

The asymptotic result of the two-point function for t — > 00 has been studied in Ref. [ll[. Calling G n — 
($(n, oo)$(0, 00)), in the case ho = 00 it has been found 



1 



G, 



2 n h™ ' 

4? cos[n arccos(Zii)] 



for/ii > 1 
for hi < 1 



(70) 



instead for ho = 



2" 



■ cosh 



(, + l)ln(i±^pi 



for hi > 1 , 
iovhi < 1 , 



(71) 
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that for large n decay exponentially with n. It has been shown that G n is decaying exponentially with n for general 
ho and hi |11|, although closed forms are not available. The exponential decay is the prediction of CFT that (maybe 
surprisingly) applies even far from the critical point hi = 1. 

The time dependence of two-point function has been studied in Ref. Q by means of exact diagonalization of the 
model with open boundary condition at the two ends r = 0, L. The results of interest for this paper are 

• The connected two-point function of a z in the thermodynamic limit and at the critical point is 

2 - ! 1 ^ J 2'-+i( 4 *) J 2-i( 4 *). (72) 

which is valid both for ho = 0,oo. Neglecting fast oscillations, (a^.(t)aQ(t)) c increases as r 2 for r < It and then 
it drops almost immediately to 0. 

• The autocorrelation function (not connected) of a z at different times is 

K(*iK(* 2 )) = Jg(2*2 ~ 2h) - i[/(t 2 + ti) ± g(t 2 - ti)} , (73) 

where f(x) — J 2 (2x) + Jq(2x), g(x) = J 2 (2x) — Jo(2x) and the sign + (— ) refers to ho = (ho = oo). 

Let us comment on these results in view of the general understanding we found so far. For t < t* all the connected 
correlation functions are zero, in agreement with CFT. All the oscillation terms of the asymptotic form cos(4£) are, 
as we discussed in the previous section, a lattice effect. Concerning the a z correlator, the r 2 dependence for t > r/2 
is a consequence of the fact that a z is not primary. The same is true for the two-time correlations function of a z that 
also decays as a power law of t — s for large times, instead of the exponential prediction by CFT for primary field. 



V. HIGHER DIMENSIONS 



Until now we just considered one-dimensional systems. Despite the fact that in low dimensions the effect of fluctu- 
ations is more pronounced making the physics highly non-trivial it is desirable to have results in higher dimensions 
as well. The method presented in Sec. [TT] to obtain the non-equilibrium dynamics of a quantum model close to 
a critical point from the critical behavior of a system confined in a slab geometry applies to generic dimension d 
through the study of the hamiltonian (T?)). Its analysis proceeds via field-theoretical RG that may provide the all 
scaling quantities of the model in an expansion close to the upper critical dimensions (u.c.d.), that is D = d + 1 = 4. 
Above the u.c.d. mean-field (or gaussian) results are exact, with logarithmic correction at the u.c.d.. For dimensions 
lower than the u.c.d., the scaling quantities are obtained as series in e = 4 — (d + 1). Thus for the time-evolution 
problem the simple mean-field solution represents an exact scaling result (a part log corrections) for the physically 
relevant three-dimensional case. An alternative method to attack analytically the hamiltonian (TjQ) is to consider an 
iV component field 4> and taking the limit N — > oo, but this will not be employed here. 



A. Dirichlet boundary conditions: the two-point function 



The Z?-dimensional slab geometry with Dirichlet boundary conditions has been the subject of several investigations. 
The two-point function has been calculated at the first order in e expansion in Ref. [431 ]. The gaussian two-point 
function, with partial Fourier transform in the parallel directions reads [43| 

1 / p -b(zi-z 2 ) i p—b(z2—zi) _ p -b(z 1 +z 2 ) _ p b(zi+z 2 ) \ 

G{p,Zi,z 2 ) = ^(e- b ^-^-e- b ^+^ + - ±£ -^—1 ? j , (74) 



with b = \/p 2 + m 2 . We are interested in the case where L = 2to, Zi — z-i — r that we will analytically continue to 
t — ► tq -\- it, and for computational simplicity we will restrict to the massless case m — 0. Thus in real space and 
imaginary time, we have (p = |p|) 



G(r,r) 



[ d3 P r -iv r 1 (, _ -2 P r , 2(l-cosh(2pr)) \ _ 
J (2tt) 3 2p \ e 4 ^» - 1 J 



(75) 
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This integral can be performed by making the sum over all the residues coming from the denominator e 4pT ° — 1. The 
calculation is rather involved, but the final result is very simple: 

1 coth(7rr/4ro)sin 2 (7rr/2To) 

G{r,T) = —. 7- r 7 ; r . (76) 

nrr cosh{nr /2t ) ~~ cos{ttt/to) 

Continuing to real time r = tq + it we obtain 

1 coth(7rr/4T )cosh 2 (^/2T ) 

G(r,t) = —, j- r —. ; r, (77) 

ttttq cosh(7rr/2ro) + cosh(7rt/ro) 

that for t, r ^> r simplifies to 

! g^rt/ro r e 7r(t-r/2)/2T U f 0Tt<r /2, , . 

G(r,t)~ -7 7^— oc < , ' , (78) 

v ; rr e ri/To+ e w/2Tii) n/ r fori>r/2. 

We recall that this result in d = 3 is exact a part log corrections. Thus the basic structure of the two-point function 
in 3D is the same as in ID, with a characteristic time t* — r/2. Using the result in Ref. [43| it is in principle possible 
to calculate the correlation functions for d < 3 in the e expansion framework. However, this requires the analytical 
continuation of complicated functions, resulting in a quite cumbersome algebra, as we shall see in the following for a 
simpler observable. 

Eq. (|74p can be used in principle to determine the gaussian behavior in any d < 3. Unfortunately the integral 
one gets is not analytically tractable for d ^ 1, 3. A possible strategy would be to perform the analytic continuation 
before of the integral and then evaluate it through a saddle-point approximation. It is straightforward to show that 
the result obtained in this way is equivalent to what we discuss in section IV Dl where we remand for the analysis of 
1 < d < 3. 



B. Dirichlet boundary conditions: a non-trivial one-point function 

In the case of Dirichlet boundary conditions, the order parameter profile in the slab geometry is trivially vanishing. 
However not all the one-point expectation values are zero. Let us consider as a typical example the operator O = 4> 2 
that has been calculated at the first order in e expansion in Ref. j44|. It has the scaling form 

(0(z,L))~L- d +V»H(z/L), (79) 
where v is the correlation length exponent. In D = 4, the function H(x) is [44| 

H(x) = ^—-^, (80) 

sin {TTXJ o 

that, using L = 2tq and continuing to z = r = tq + it, leads to the real time evolution 

(Off)) - (Oft = 00)) ~ T - 2 JL r -2 -rt/iha / 81 j 

V K " V V ;/ cosh 2 (7rf/2T ) V ' 

i.e. (0(t)) approaches its asymptotic value exponentially, as it as been found for all primary operators in d = 1. 
The one-loop result for the O(N) model in D = 4 — e dimensions is [44j 



H(x) 



sin (nx) 3 



( 1 + In - e[C'(2, x) + ?(2, l-x)]+ const , (82) 



where const stands for terms do not depend on x and C'( a ? x ) = d a £(a, x) and ^ is the generalized Riemann function. 
This example clearly shows that within the e expansion, the scaling functions contain logarithmic contributions that 
originate from the expansion of power law as e.g. sin e irx = 1 + elnsin7ra;. To have a function with good analytical 
structure to perform the real time continuation, it is desirable to "exponentiate" such logarithms. In Ref. [44j | the 
exponentiation procedure leaded to 



2—1/1/ 

( — — " ) ±lD-2..r) + ^D-2.l-.r ) -2^D 



2 / \ D-2-l/v 

^(2-1/1/) -1^ , (83) 

6 \sm(7rx) / 
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with v = 1/2+ (N + 2)/(N + 8)e/4. Performing the analytical continuation we have the sum of two exponentials, 
and in any dimension d < 3 the second one has a largest "relaxation time" that hence is dominating for large t. 
Considering only the second term we have 

(0(t > r )) - (0(t = oc)) oc e -(d-i-i/")Tt/2T _ (84) 

Note that the subleading term e -( 2 - 1 / l/ ) 7rt / 2T o j g multiplied by a \ogt term arising from the £ function. 

This example put forward the idea that (at least for Dirichlet boundary conditions) the exponential relaxation of 
the one-point functions is not only a property of one-dimensional systems, but holds in any dimension (with eventually 
log corrections) with relaxation times related to the scaling dimensions of the operator. 



C. Fixed boundary conditions: the one-point function 



The case of the extraordinary transition, that corresponds to fixed boundary conditions, has been considered in 
Ref. (45[. The magnetization profile can be written as [45| 

,. . 2K i dn(2K i z/L) 

<K*)= T ; * ~TT\ ' 85 

L sn(2iV4z/_L) 

where Ki = K(k^), where K(k) is the elliptic integral, ki the elliptic modulus that in terms of the parameter of the 
model is cj> 2 L 2 = (2_ftT 4 ) 2 (2fc| — 1), and sn(x) and dn(x) and the Jacobi functions. Continuing to real time, and using 
the properties of the Jacobi functions we obtain 

(f>(t) = y/i - k 4 cn(K 4 t/T , 1 - fc 4 ) , (86) 

that, contrarily to all the other cases we have considered, is oscillating and not exponentially decreasing. 

To our knowledge there are no result in the e-expansion for the magnetization profile. However, using "local- 
functional methods" [46| it has been obtained an approximate profile in D = 3 that involves, as in mean-field theory, 
Jacobi elliptic functions. It can be easily continued to real time via z — > tq + it and again one finds an oscillating 
behavior with time. The method exploited in Ref. 46] can be used in any dimensions, obtaining an always an 
oscillating behavior with a period that diverges as D approaches 2, recovering Eq. ([7|). All these calculations are 
essentially mean field and we do not know how the inclusion of fluctuations changes them. It can possible that for the 
extraordinary transition (i.e. fixed initial conditions), the exponential decay founds at D — 2 is more an exception 
rather than a rule, because of the simple analytic structure of the trigonometric functions in the complex plane. 
Another possibility is that fluctuations destroy these oscillations. Only a complete analytical calculation (e.g. in large 
N) can help in understanding this point. 



D. A real-time solvable model 



As for the one-dimensional case, it is worth to check the results coming from the analytical continuation of large 
imaginary-time with exactly solvable models. The simplest (and probably one of the few) model solvable in generic 
dimension is the generalization of the hamiltonian (|38|) to a d-dimcnsional hypercubic lattice. The solution of such 
model proceeds via Fourier transform as in one dimension. The final result is simply give by Eq. (|53p with the 
replacement dk — > d d k, i.e. 

(<p r (t)M*)) ~ (<Pr(P)M0)) = / 7^ e — T727i ^—H-li. (87) 

First of all let us note that this expression in virtually identical to Eq. (|75p if we take r oc m^ 1 — > 0. In fact, taking 
To — ► in Eq. (|75|) . only the third term matters, since it is O^q 1 ) relative to the first two. Taking z = r + it and 
b — fl p we get G(p) oc t^ 1 J S!p 2 (l — cos2fl p t) that is Eq. (|87|) for m — > 00. 

To understand the general features let us consider in details the conformal evolution (f2/- = w|k|) from a disordered 
state (Slofc = wo). The derivative of the two-point function is (k — |k|) 



d t ((p r (t)(p (t)) = 2m Im 



^dj, gi(k-r— 2kvt) 

(2n) d k 
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Except for d = 1, this can be done analytically only in d = 3, where we can write it as 

kdk j d9 sin 9e lk( - rcos 9 ~ f dk(sm(kr)/r)e 2ikvt ~ (l/r)d(trt - r/2) . (89) 



Integrating with respect to t we get zero for t < r/2v and 1/r for < > r/2v. 
For general d we have 



d t (cp r (t)cp (t)) cx2™ Im / k d - 2 dk / d0(sin0) <i - 2 e i * i(rcoB ' 9 - 2, ' t) . (90) 



By a saddle point argument, we can assume that the dominant behavior comes from 9 close to zero, so the 9 integral 
gives 



leading to 



d-2 e - lk re-/2 „ ( fer )-(d-l)/2 ; ( 91) 



dk— J^-l—e^- 2 ^ ~ r-( d - 1 )/ 2 (2«t - r)-( rf - x )/ 2 e(^ - r/2) . (92) 
{kr)( >i 



Integrating with respect to t, we get zero for t < r/2v, as expected by causality, and 

r -(d-i)/2( 2ut _ r j( 3 -«0/2 ) (93) 

for t > r/2v. It is interesting that this gaussian correlation function blows up at large t only for d < 3, when we 
expect the fluctuations to become important. It would be interesting to study this in the (j) A theory for large N, by 
replacing 4> 4 by 3{(f> 2 )(f> 2 in the usual way, where now (4> 2 ) depends on t and is calculated self-consistently. 



VI. PHYSICAL INTERPRETATION AND DISCUSSION 



In this paper we studied in general the non-equilibrium unitary dynamics that follows a sudden quantum quench. 
We showed that if the hamiltonian H governing the time evolution is at a critical point, while Ho (i.e. the one for 
t < 0) is not, the expectation value of a class of operators (primary ones in CFT) relaxes to the ground-state value 
exponentially in time with universal ratio of decaying constants. We also found that connected two-point functions 
of operators at distance r are vanishing for t < r/2v, while for t > r/2v reach exponentially fast a value that depends 
exponentially on the separation, in contrast with the power laws typical of equilibrium configuration. 

We also considered the real-time dynamics of simple exactly solvable models and we found that several of the typical 
characteristics of the critical points still hold. Roughly speaking, critical points are not special as far as quenching 
dynamics is concerned. In fact, also for gapped systems, connected correlation functions vanish (or are strongly 
suppressed) for t < r/2v and for asymptotic large times resemble those at finite temperature despite the fact that the 
whole system is in a pure state. Several other examples in the recent literature (see e.g. 0, [H, 0, [H, H^, |4l| ) gives 
further evidence that these two effects are actually true in general, at least in the realm of exactly solvable models 
considered so far. 

In the following we give a simple interpretation of these two features separately trying to understanding their 
physical origin. 



A. The horizon effect 



The qualitative, and many of the quantitative, features found for the time evolution of correlation functions may 
be understood physically on the basis of a picture we first introduced in Ref. [39( to describe the time evolution of the 
entanglement entropy. Later we generalized it to correlation functions in Ref. [33| and it has been largely adopted 
thereafter dill US S3- 

We emphasize that such scheme is not an ab initio calculation but rather a simplified picture which allows us to 
explain physically our findings. The initial state \ipo) has an (extensively) high energy relative to the ground state 
of the hamiltonian H which governs the subsequent time evolution, and therefore acts as a source of quasiparticle 
excitations. Those quasi-particles originating from closely separated points (roughly within the correlation length 
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£o = tRq of the ground state of Hq) are quantum entangled and particles emitted from far different points are 
incoherent. If the quasiparticle dispersion relation is E = flk, the classical velocity is v& = Vfefife. We assume that 
there is a maximum allowed speed v m — max^ |vk|. A quasiparticle of momentum fc produced at r is therefore at 
r + Vfei at time t, ignoring scattering effects. This is the only physical assumption of the argument. Scattering effects 
are not present in the theories considered so far, but as evident from the argument outlined below they can play a 
role for only t > r/2v m (allowing, perhaps, for a renormalization of v m by the interactions). 

These free quasi-particles have two distinct effects. Firstly, incoherent quasi-particles arriving a given point r from 
well-separated sources cause relaxation of (most) local observables at r towards their ground state expectation values. 
(An exception is the local energy density which of course is conserved.) Secondly, entangled quasi-particles arriving 
at the same time t at points with separation |r| 3> £o induce quantum correlations between local observables. In the 
case where they travel at a unique speed v (as in CFT), therefore, there is a sharp "horizon" or light-cone effect: the 
connected correlations do not change from their initial values until time t ~ |r|/2u. In the CFT case this horizon effect 
is rounded off in a (calculable) manner over the region t — \t\/2v ~ To, since quasi-particles remain entangled over this 
distance scale. After this they rapidly saturate to time-independent values. For large separations (but still <C 2vt), 
these decay exponentially ~ exp(— 7ra;|r|/2ro). Thus, while the generic one-point functions relax to their ground-state 
values (we recall in CFT this relaxation is exponential ~ exp(— nxvt/To)), the correlation functions do not, because, 
at quantum criticality, these would have a power law dependence. Of course, this is to be expected since the mean 
energy is much higher than that of the ground state, and it does not relax. 

This simple argument also explains why for the case of a semi-infinite chain the relevant time scale is r/v rather 
than r/2v, since one of the two particles arriving in r has been reflected from the end of the chain. This has been 
also stressed in Ref. [4l|, in the study of the time dependence of entanglement entropy of finite chains with open 
boundary conditions. 

All our results are consistent with this picture as long as the quasi-particles are assumed to all propagate at the 
same speed, resulting from a "conformal" dispersion relation f2/- = v\k\. However, it is very simple to generalize this 
picture to different dispersion relations, taking into account that each particle propagates at group velocity Vk = in- 
appropriate to the wave number k. In this case the horizon effect first occurs at time t ~ |r|/2i> m , where v m is 
the maximum group velocity. If v m occurs at a non-zero wave number, it gives rise to spatial oscillations in the 
correlation function. Thus again connected correlation function are expected to be strongly suppressed for t < t* 
and start developing only after t*. In the case of a general dispersion relation we do not have a proof, beyond the 
stationary phase approximation, that the connected correlation functions remain constant up to this time, but one 
would expect it on the grounds of causality. (The proof in the case of a relativistic dispersion relation uses Lorentz 
invariance in an essential way.) However, because there are also quasi-particles moving at speeds less than v m , the 
approach to the asymptotic behavior at late times is less abrupt. In fact, for a lattice dispersion relation where Q' k 
vanishes at the zone boundary, the approach to the limit is slow, as an inverse power of t. A similar result applies to 
the 1-point functions. This is consistent with the exact results obtained here and elsewhere. 



B. The large time limit and the generalized Gibbs ensemble 



The existence and the understanding of the asymptotic state resulting from the evolution from an arbitrary state is 
one of the most-interesting problem in statistical mechanics. A robust theory able to predict this state ab-initio still 
does not exist. A currently popular idea is that for late times the system (or rather macroscopically large subsystems) 
'look like' they are in a thermal state, despite the fact that the actual state of the whole system is pure. A common 
belief is that a region of dimension r can be thermalized by the infinitely large rest of the system which acts as a bath 
(see e.g. 0, IH, H?})- But this intriguing idea is not sufficient to give the value of the resulting effective temperature. 

A major step toward the clarification of the properties of the asymptotic state has been made by Rigol et al. [ll|. 
In fact, it was conjectured that if the asymptotic stationary state exists, it is given by a generalized Gibbs ensemble 
obtained by maximizing the entropy S = — Trplogp, subject to all the constraints imposed by the dynamics [l5| . 
Consequently, denoting with J m a maximal set of commuting and linearly independent integrals of motion, the density 
matrix is 

p = z -i e -j: m i m i m) z = Tre -J2 m i m i m _ (94) 

We note that such a density matrix describes a pure state only if the model under consideration is integrable, i.e. if 
the number of integral of motions equals the number degrees of freedom. If there are not enough integrals of motion, 
p corresponds to a mixed state and it is not clear to us to which extent it can describe the pure state resulting from 
the time-evolution. The values of the Lagrange multipliers l m are fixed by the initial conditions: 



Tr I m p = {I m ) t =o 



(95) 
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In the following we will take this generalized Gibbs ensemble as a postulate and we will show how it nicely and 
naturally explains the "effective temperature" effect observed for large times. However we stress that there is still no 
proof for this assumption that, to our opinion, cannot be considered on the same fundamental level as the thermal 
Gibbs ensemble. 

Let us consider the chain of harmonic oscillators of the previous section as a typical example. We will soon see that 
most of the features are quite general. In this case the natural choice for an infinite set of integral of motion is the 
number of particles with momentum k, i.e. n k = A k A k . Most other observables can be written in terms of these, i.e. 
H = ^>2 k Qferifc. Consequently the expectation value is given by 

(O) t=00 = Ti Op = Tr OZ- 1 e- ^" lknk , (96) 

that can be seen as a thermal density matrix with a k dependent effective temperature given by 

/3 cff (fc)fi fc = l k . (97) 

Thus an effective temperature already appeared. Note that this state can still be pure because such a temperature is 
k dependent. 

To fix lk we need (n k )t=o- From Eq. (|50|) we get 

w„-<4^-4-i(£ + £)-i. (98) 

The calculation then proceeds as for a thermal distribution 

d °° 1 

(n k ) p = Trn k p = - — \nZ , with Z = Tr e~ E * l " n " =T[z2 = IT 1 _ e -i k ' (") 

so that 



k n k =0 



Wp4EMi-^) = ^. (wo) 

k 



From this e lk = 1 + d k 2 and 

n k 

At finite temperature the correlation function in momentum space is 



/M*) = . (101) 



9 9 / \ 9 1 _l p—0Qk 

(fPkf-kh = {^{A k + Ai){A_ k + Al k )) p = — ((AlA_ k ) p + (A k Al k ) p ) = — 1 _ ^ , (102) 



L k ii k 

that substituting the previous result for (3 c s gives 



p{k) t =oo = (<Pk<P-k)t=oo = 7^(1 + 24) = "L + "° fc , (103) 

that is exactly the time-average of Eq. (|54[) reproducing, after Fourier transforming, the well defined correlation in 
real space for t — > oo. 

So this generalized Gibbs ensemble correctly reproduces the exact diagonalization of the model and gives also 
few insights more. In fact, in the limit m,Q — > oo the effective temperature is independent from k and m obtaining 
/3 e fi = 4/mo, explaining a posteriori the simplicity of the results in this case. Note instead that for arbitrary nig and 
m = 0, i.e. conformal evolution, f3 e s(k) is a function of fc. In this case, the large distance properties of correlation 
functions are described by the mode with k = 0, for which, independently of m we get f3 c a{k = 0, m = 0) = 4/mo, 
consistently with the previous findings and general expectations. Furthermore we can conclude that the large r 
asymptotic behavior is always governed by the effective temperature /3 e g(k = 0) = 2(log(|m — m|/(m + m))/m. 
Another interesting feature is that /3 e s(k = 0) gives the asymptotic behavior of the correlation functions of all those 
observables that are effectively coupled with the zero-mode. In the opposite case the relevant temperature is the 
largest /3 e s(k) with the k mode coupled with the observable. Another interesting feature is that on a very rough basis 
one can be tempted to assume that that /3 e s(k) is directly related to the excess of energy of the mode k generated by 
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the quench, but this is not the case. In fact (H)t=o — J2k^kdk that is different from /3 eff , being the same only in the 
limit d^ 2 — > 0, i.e. mo — ^ oo. 

Clearly the same reasoning of before applies every time we consider a model that can be diagonalized in momentum 
space with a proper choice of the quasi-particles, i.e. every time that H = Q,kA k Ak for some Af.- This means that 
the excitations are non-interacting. As far as we are aware all the applications of the generalized Gibbs ensemble to 
the date only concern this kind of models [HI, Hi| , but there are few numerical hints suggesting that can be true more 
generically [291 ]. 

Now we outline how we imagine the generalized Gibbs ensemble given by Eq. (|94p could be used to justify the 
effective temperature scenario for large time for any integrable system, i.e. with a complete set of integrals of motion. 
The hamiltonian can be written as H — Y] m o, m I m , with some a m eventually zero. Thus one can think to an in- 
dependent effective temperature /3 e fj(m) = / m /a m , but this temperature does not give directly the behavior of the 
correlation function for large distance because in general the integral of motions are not diagonal in fc-space. Thus we 
can only conjecture that the correlation functions of a given operator O(r) are governed by the largest p c g(m) with 
to among the integral of motions to which O(r) effectively couples. We stress that this is only a crude argument and 
we are still not able to put it on a firmer basis. 



C. Discussion and open questions 



We presented a quite complete picture on the time evolution of a quantum system after a sudden quench of one 
hamiltonian parameter. Despite the fact that a lot of work has been done, still more is left for future investigation. 

The first problem that must be addressed is the real-time dynamics of effectively interacting systems. In this 
direction Bethe Ansatz solvable models like the Lieb-Liniger gas and Heisenberg spin chains are among the best 
candidates for an analytical approach. This would clarify how a non-trivial two-body scattering matrix can modify the 
time-evolution (only inside the light-cone) of the correlation functions and whether the crude argument we outlined 
for an "effective temperature" for the long-time state is valid. To clarify this point also the study of (boundary) 
integrable massive field theories can be of some help. 

Another approach, currently under investigation: 49], is the direct analysis of the perturbative expansion for the 
correlation functions in a A</> field theory. This has some simplifying features in the large too limit which may allow 
it to be resummed to all orders. It is very important to get non-trivial results for such a non-integrable interacting 
theory. 

Numerical computations for non-integrable systems should also be performed to understand if and eventually to 
which extent the generalized Gibbs ensemble picture is valid beyond integrability. Some numerics concerning this 
point are already available [22|, HH, Hil but still a clear scenario is not emerged. Among non-integrable models special 
care must be given to disordered systems because of the non existence of a speed of sound as a consequence of 
Anderson localization 49]. Some insights can be obtained from those models whose equilibrium behavior can be 
analytically obtained by means of the strong-disorder renormalization-group Even in this case time-dependent 
DMRG can help a general understanding. The time-evolution of the entanglement entropy 4l| already revealed that 
in these systems there is no sign of the light-cone (as easily predictable because of the absence of speed of sound) and 
the effective motion of "quasi-particles" is more diffusive rather than ballistic. Furthermore the results at the largest 
available times does not look thermal at all, but this can be also due to the time window accessible with numerics. 

Another very interesting question is what happens at finite temperature and under what conditions a system can 
equilibrate. It should be relatively simple to generalize the results for all the "quasi-free" models already considered 
at T = 0. Work in this direction is in progress, but it will be almost impossible to consider the same problem for 
more complicated models. 

Finally it is also important to understand the role played by the initial state. We assumed always a translational 
invariant one with short range correlations. Thus one natural modification consists in taking a state that is only 
locally different from the actual ground-state. This problem is known in the literature as a "local quench" and it has 
already been considered to some extent (see e.g. [FJ, [EU, [53[), but still a general picture as the one outlined here 
for glob al quenches does not exist (this is also interesting for eventual connections with quantum impurity problems, 
see [54| as a review). Another natural modification of the initial state is one with long-range correlations, such as a 
critical one and let it evolve with another critical hamiltonian. This has been done for the Luttinger liquid [261 ] and 
the results show the typical functional dependence of a light-cone scenario (i.e. everything depends only onil 2vt) 
but the long-time correlations decays as power-laws, with exponents that are different from equilibrium ones and can 
be predicted by the generalized Gibbs ensemble. 
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